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In this talk, we discuss some general properties of particle production 
in a field theory coupled to strong time dependent sources, and techniques 
. to compute the spectrum of the produced particles in such theories. We 

*■■»«. ' also discuss the application of these results to the description of hadron 

' or heavy ion collisions in the Color Glass Condensate framework. 
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1 Introduction 



At high energy, all the internal timescales of a hadron are time dilated. There- 
fore, more and more soft fluctuations - carrying a smaller and smaller fraction 
y\ ' x of the hadron momentum - become long-lived and become relevant in interac- 

tions with another hadron. On the contrary, on the timescales relevant for such 
an interaction process, the large x partons can be seen as completely frozen de- 
grees of freedom, whose only role is to act as sources that radiate more small x 
gluons. Moreover, the small x modes will eventually have an occupation number 
larger than unity, and undergo recombinations - a process known as saturation 

m 

In the Color Glass Condensate (CGC) framework [313], one thus divides the 
degrees of freedom of a hadron in static color sources - described by a density 
denoted p - that represent the large x partons, and dynamical gauge fields that 
represent the small x partons. The CGC can thus be seen as an effective theory 
of gauge fields coupled to external sources. The details of this separation of 
degrees of freedom can change with the separation scale, but this should not 
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affect physicaf quantities. This leads to a renormalization group equation - 
known as the JIMWLK equation [3] - that governs the evolution with x of the 
distribution W[p] of hard color sources. 

In these proceedings, we consider the collision at high energy of two hadrons 
(or heavy ions) described in the CGC framework. We assume that the distribu- 
tions of hard color sources that describe the two projectiles are known, and we 
address the question of calculating physical observables in given configurations 
of the two sources. Moreover, we will consider only the regime where the two 
projectiles are saturated, in which the two sources are strong - both of 
order g~ 2 . 

As stated before, one must consider an effective theory described by the 
following Lagrangian, 

£ = -\^ V v + \ (Ji + > (i) 

with currents given at lowest order in the sources by 

J? = g5» + 5(x-) Pl (x ± ) , J£ = gS»-8(x + ) Pi (x x ). (2) 

These currents must be covariantly conserved - [D^, Jf 2 ] = - which leads to a 
feedback of the gauge field on the currents. Therefore, in general, the eqs. (J2J) get 
modified by corrections of higher order in pi^- An important observation is that 
the strength of the sources lead to non-perturbative effects, in the sense that 
an infinite set of diagrams must be summed in order to calculate a quantity 
at a fixed order in g. However, the large strength of the sources has also a 
valuable consequence: as we shall see later, the leading order is dominated by 
tree diagrams only, and it can be studied by classical methods. 



2 AGK cancellations 

Following U, let us first consider the theory of a real scalar field coupled to 
strong sources. Most of the structural properties we want to discuss can indeed 
already be studied in this simpler framework. Some general results - that are 
to a large extent independent of the details of the theory under consideration 
can be obtained by considering the generating function for the probabilities P n 
of producing a given number n of particles : 

00 

T[z) = Y i P n z n . (3) 

In 0], we have proven that, if we denote b r /g 2 the sum of all the cut connected 
vacuum- vacuum diagrams 1 with exactly r cut lines, the logarithm of J-{z) reads 

oc 

hf(z) = -j^6 r (/-l). (4) 

" r=l 

lr The explicit factor l/g 2 represents the natural order of these cut diagrams when the 
sources are strong. 
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From eqs. (0 and one gets the following expression 2 for the probability P n : 
Pn = e-^±' ^ (5 ) 

p=0 F nH hr p =n y 

In this formula, the index p is the number of connected cut subdiagrams from 
which the n particles are produced: r\ particles are produced by the first of 
these p subdiagrams, ri by the second one, etc... If we further expand the 
prefactor exp(— J2 r b r /g 2 ), we obtain : 

<3=0 r p=0 riH hr p =n 3 

In this expression for P n , the index q is the number of connected subdiagrams 
that are not cut (this can be seen from the fact that it comes from the absorptive 
correction, whose sole role is to preserve unitarity). 

The term of fixed p and q in this formula can therefore be interpreted as 
the probability of producing n particles from p + q connected subdiagrams, p 
of which are cut and q of which are not cut. By summing this probability 
over n, one "integrates out" some degrees of freedom in order to keep only the 
information about the probability of having p cut subdiagrams and q that are 
not cut, regardless of the number of produced particles. This probability reads : 

Summing this expression over q from to oo, one finally obtains the probability 
of having p cut subdiagrams : 

One therefore sees that the number of cut subdiagrams has a Poissonian distri- 
bution, with an average of (n) t — Yl r ^r/g 2 - Eqs. (7J and (JSJ) are the essence 
of the Abramovsky-Gribov-Kancheli cancellations UJ. As one can see from the 
above derivation, they are simply a consequence of the factorization of a generic 
diagram in terms of its connected subdiagrams. Therefore, we expect them to 
be much more general than the context of reggeons field theories in which they 
have first been discussed. Another point should also be obvious at this point : 
in order to obtain the eqs. Q and iJSJ, one has "integrated out" the number n 
of produced particles. By doing this, the infinite sequence 61,62,63, • • • has re- 
duced to the single combination ^ r b r . This means that a lot of the dynamical 
information about the theory under consideration has been lost in this process, 
and that there are certain questions that cannot be answered anymore by the 
sole knowledge of the lZ p . q 's. 

2 A model for the coefficients b r has recently been proposed in 
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3 Generating function 



Let us now consider the generating function T(z) per se, and discuss what it 
would take to calculate it. The following results on this question have been 
established in j3] : 

(i) Diagrammatically, T(z) is the sum of all the cut vacuum-vacuum dia- 
grams, where each cut propagator is weighted by a factor z. 

(ii) At leading order, the derivative of \a.T{z), T' (z) j T{z) can be calculated 
from a pair of solutions 4>±(z|x) of the classical equation of motion. For 
a scalar field theory with a cubic coupling and a source j, the classical 
EOM reads 

(□, + m 2 )$±(z\x) + |$|0|z) = j(x) . (9) 

The expression of T\z)jT(z) is simpler if written in terms of the Fourier 
modes f±~\z\x°,p) and f±(z\x°,p) of these classical fields, 

*±(z\x) = I {2 ^ 3 P 2Ep {.f { ± +) (z\x ,p)e-^ + fi- ) (z\x ,p)e + ^} , 



(10) 



and reads 
T\z) 



Hz) 



(iii) The two solutions $± of the classical EOM must obey the following bound- 
ary conditions : 

fi +) (z\x° = -oo, p) = , f ( S\z\x° = -oo, p) = , 

f M (z\x° - +co,p) - zfi +) (z\x° = +oo,p) , 

f { + ~\z\x° = +oo,p) = zf { T ) {z\x Q = +oo,p) . (12) 

(iv) From unitarity, J-(l) = ^2 n P n = 1. This property serves as the initial 
condition for going from T\z)jT{z) to T{£), by writing : 

^(z) = exp||dr^l . (13) 



Unfortunately, finding the pair of solutions of the classical equation of motion 
that obey the boundary conditions of eqs. I|12H is a difficult numerical problem, 
that has not yet been studied in this context. However, as we shall see in the 
next section, one can obtain from here a formula for the average multiplicity 
which is much easier to evaluate numerically. 
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4 Average multiplicity at leading order 



4.1 General method 

The moments of the distribution of multiplicities, in particular the average mul- 
tiplicity (n) = ^2 n nP n , enjoy a special status because extra simplifications 
occur in their calculation. Let us consider the multiplicity since it is the sim- 
plest one. One can get the multiplicity from the generating function as 

(n) = = 5$ ■ (14) 

(The second equality is of course simply due to = 1.) Therefore, calculating 
(n) is a special case - with z = 1 - of the calculation of T'(z)jT[z). At z = 1, 
the third and fourth of the boundary conditions in eqs. <|12[) simply reduce to 

= +oo,p) = fi + \l\x° = +oo,p) , 
f { + -\l\x° = +co,p) = f ( _-\l\x° = +oo, p) , (15) 

which means that the fields $+ and $_ are equal at large positive times, as 
well as their first time derivative. Since the classical equation of motion is 
deterministic, this implies that these two classical fields are equal at all times. 
The first two boundary conditions in eqs. (|12|) then imply that <£>± are vanishing 
at large negative times, as well as their first time derivative. 

Therefore, at leading order, the multiplicity (n) is given by eqs. (|1(J|I and 
with <i> + = the retarded solution of the classical EOM with a null initial 
condition at xq = — oo. One should emphasize the following: it was obvious 
from the beginning that the multiplicity at leading order would be somehow 
related to solutions of the classical equation of motion - since it involves only 
tree diagrams at this order - but it is a non-trivial result that this is with 
retarded boundary conditions. The retarded nature of the boundary condition 
is crucial in practice, in order to solve this problem numerically. 



4.2 Gluon multiplicity at leading order 

Going back to QCD, it is straightforward to generalize the previous results to 
the case of gluon production. The inclusive gluon spectrum is given at leading 
order by 



d(n g iuons) i f 

e p d3p — = Y,J xy eip<x ' v) n ^y ^ ■ A *w ^ ■ A ^y) 



(16) 



In this formula, A 11 (x) is the retarded solution of the classical Yang-Mills equa- 
tions - in the presence of the sources px t 2 - and with A^t = and d°A 1 ^ = at 
xq = — oo. This problem was solved numerically in The diagrams 3 that are 

3 The propagators that appear in this diagrammatic representation are not Feynman prop- 
agators, but retarded propagators. 
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Figure 1: Left : space-time representation of the diagrams involved in the cal- 
culation of the gluon multiplicity at leading order via eq. I|16l) (one factor of 
A 1 ^ is represented in black and the other in grey). Right : the resulting gluon 
spectrum - the plotted quantity is dN/d 2 k±. 

resummed by eq. (|ltj|) arc represented in figure ^ as wei l as the resulting gluon 
spectrum. 

4.3 Quark production at leading order 

Similarly, the production of quarks at leading order has been considered in 
[B]. Note that we use "leading order" somewhat abusively here, since quark 
production is strictly speaking a Next-to-Leading Order effect - indeed, in the 
saturated regime, the number of produced gluons scales like g~ 2 while the num- 
ber of produced quarks scales like g°. For quark production at leading order, 
one must use the following formula, 

E P d3p /lo = ^3 J J e*-<*-»> (i0 x -m)(i0 v +m) %(x)Mv) > (17) 

where ip q {x) is the retarded solution of the Dirac equation - with the classical 
field obtained in the calculation of gluon production in the background - with 
a free negative energy spinor, v(q)e lq ' x , as the initial condition. The diagrams 
involved in eq. Ijl7|l are sketched in the left of figure Also represented in 
the right part of this figure is the resulting quark spectrum, for various quark 
masses. 
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Figure 2: Left : space-time representation of the diagrams involved in the cal- 
culation of the quark multiplicity at leading order. Right : the resulting quark 
spectra for various quark masses. 

5 Further developments 
5.1 Gluon multiplicity at NLO 

In fact, the production of quarks is one of the pieces - the simplest one - that 
contribute to particle production at NLO, i.e. at order g°. In ^Q, we have 
detailed the principles of a full NLO calculation of the particle yield, in the 
case of scalar fields. Of course, things will be more complicated in QCD with 
gluons, but one crucial property of the result will survive : the calculation of 
particle production at NLO can be done from retarded solutions of the classical 
EOM and retarded solutions of the EOM of a small fluctuation on top of the 
classical held. The crucial point again is the fact that these objects are needed 
with retarded boundary conditions, which means that it is a problem which is 
tractable numerically in a straightforward way. 

Two types of topologies, sketched in figure [31 contribute to the gluon multi- 
plicity at NLO 9 . One of them (left diagram) is very similar to that already 
encountered in quark production - it corresponds to the production of pairs 
of gluons, and involves retarded solutions for the equation of motion of small 
gluonic fluctuations on top of the classical field. The diagram on the right can 
be seen as a 1-loop virtual correction to the classical field - the field in the 
complex conjugate amplitude remaining the tree-level one. It was shown in |3] 
that the latter contribution can also be expressed in terms of retarded solutions 
of the small fluctuations equation. Therefore, the result according to which the 
inclusive multiplicity can be calculated from retarded solutions of some equa- 
tions of motion remain true at NLO. To this diagram with a virtual gluon loop, 
one must add two similar diagrams, respectively with a quark loop and a ghost 
loop (only in gauges that have ghosts for the latter). 

One additional issue arises when one considers these one-loop corrections 
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Figure 3: Space-time representation of the diagrams involved in the gluon mul- 
tiplicity at NLO. 

to the gluon yield: some of the contributions have a divergence of the form 
a s J dx/x which is reminiscent of the divergences already resummed by the 
JIMWLK evolution of the distribution of sources W^pi] and W\p2\ for the two 
projectiles. For the CGC framework to be self-consistent, one must prove that 
these divergences that appear in the gluon yield for fixed p\ and pi can be 
absorbed in the evolution of the W^pi^] [HI- 

5.2 Survival probabilities 

One can also consider exclusive processes in this framework. For instance, in- 
stead of the plain - fully inclusive - probability P n of producing n particles, one 
may define a probability P n (Sl) of producing n particles in a certain region Q 
of the phase-space and none outside of J7, and construct a generating function 
for these exclusive probabilities, Fn{z) = X^^Lo Pn(ty z n ■ O ne can show JHi 
that this generating function can in principle be calculated at leading order by 
methods that are similar to those described in section[3| modulo two differences. 

(i) The boundary conditions for the two classical fields in terms of which on 
can express Tq(z) / ^(z) read : 





l*° = 


-oo, p) 


= 0, f ( -\z\x° = - 


•oo, p) = , 


f M (z 


\x° = 


+oo,p) 


= z n(p)f ( h +) (z\x° = 


+oo,p) , 




\x° = 


+oo,p) 


= z£l(p) f ( r\z\x° = 


+oo,p) , 



where f2(p) is a function which is 1 in the region f2 and zero outside. 

(ii) The "integration constant" ^h(l) is n o longer unity. In fact, this quantity 
is the total probability of not producing particles outside of SI. It will 
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appear as a prefactor in all the probabilities P n (f2), and it can therefore 
be interpreted as a survival probability for the empty region outside fi. 



6 Conclusions 

Multiparticle production in field theories coupled to external time-dependent 
sources - e.g. in the Color Glass Condensate framework for hadronic collisions 
at high energy - exhibits some non-trivial features when these sources are as 
strong as the inverse coupling: even in the weak coupling regime, one must 
resum at each order an infinite set of diagrams. Quite generically, one recovers 
in this type of model the combinatoric relations among the probabilities of cut 
subdiagrams that lead to the AGK cancellations. 

At leading order, both the generating function of the probability distribution 
and the average multiplicity can be calculated from solutions of the classical 
equation of motion in the presence of the external sources. However, while 
these solutions must be found with complicated boundary conditions in the 
case of the generating function, the problem can be reduced to finding retarded 
solutions in the case of the multiplicity, leading to straightforward algorithms for 
calculating the multiplicity numerically. This has been done at leading order in 
the CGC framework for the production of gluons and quarks in nucleus-nucleus 
collisions. 

The present study can be extended in several directions. One of these ex- 
tensions is the production of gluons at NLO (NLO particle production has in 
fact been studied in g], but the techniques developed there must be extended 
to gauge theories). Another extension is the study of exclusive reactions, where 
one enforces some constraints on the phase-space of the produced particles. 
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